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In this paper we discuss the incorporation of gauge links in hadronic matrix elements that describe 
the soft hadronic physics in high energy scattering processes. In this description the matrix elements 
appear in soft correlators and they contain non-local combinations of quark and gluon fields. In our 
description we go beyond the collinear approach in which case also the dependence on transverse 
momenta of partons is taken into consideration. The non-locality in the transverse direction leads 
to a complex gauge link structure for the full process, in which color is entangled, even at tree-level. 
We show that at tree- level in a 1-parton unintegrated (1PU) situation, in which only the transverse 
momentum of one of the initial state hadrons is relevant, one can get a factorized expression involving 
transverse momentum dependent (TMD) distribution functions. We point out problems at the level 
of two initial state hadrons, even for relatively simple processes such as Drell-Yan scattering. 
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I. INTRODUCTION 

*^ In the description of high energy processes involving hadrons one aims at isolating the underlying hard process in 

terms of partons, quarks and gluons. The purpose of this can be twofold. The aim might be to study the quark and 
O ,. gluon structure of hadrons, or it might be to account for the soft hadronic physics to study unknown details in the 
hard process, for instance involving physics beyond the Standard Model. In both cases one identifies a number of soft 
functions among them distribution functions f H ~^ l {x) and fragmentation functions D h ^ l (z) for quarks (i = q, where q 
is an (anti-)quark flavor) and gluons (i — g), which have a natural interpretation as the probability of finding a parton 
i with momentum fractions i in a hadron H or as measure of the number of hadrons h with momentum fractions z 
J> . in the 'decay' of a parton i. We will refine the definitions of these fractions below. To some level of accuracy, one can 
express observables such as cross sections and asymmetries in terms of these distribution and fragmentation functions. 
Going beyond the collinear treatment, one includes the dependence on transverse momenta. These momenta can serve 
as degrees of freedom in the connection between hadrons and partons very much like spin degrees of freedom. In 
the case of fragmentation, transverse momentum is in essence just the mismatch between parton momentum k and 
hadron momentum Kh, or better between the fraction of the parton momentum zk and Kh- This, at least intuitively, 
corresponds for jet fragmentation to identifying the parton momentum with an appropriately defined jet direction. 
Also for initial state hadrons one can include dependence on transverse momentum, which is the mismatch between an 
• • ' appropriate fraction of the hadron momentum xP and the parton momentum p. In a high energy scattering process, 
one is able to use the presence of a hard scale to identify parton momenta (integration variables) with accessible 
combinations of external momenta. This is best known for the momentum fractions, but it is also possible for the 
^ , transverse momenta. 

In the present paper, we remain at tree- level (to be made more explicit), which implies that the intuitive language can 
in most cases also be used at the level of the matrix elements of quark and gluon fields that constitute correlators 
which in turn are parametrized in terms of the beforementioned soft functions connecting partons and hadrons 
A complication that arises already at tree-level, is that the appropriate combinations of quark and gluon fields in the 
correlators need to be gauge-invariant combinations. This would trivially be the case for local products of fields, but 
already in the collinear treatment which studies the dependence on the momentum fractions x and z, the partonic 
field combinations are non-local along a light-like direction conjugate to the parton momentum. Then one finds that 
gluon fields with polarizations along the momentum direction, which naturally appear in a twist analysis of leading 
operators, need to be resummed to make up the required path ordered exponentials or Wilson lines connecting the 
non-local parton fields |3|. Although this involves an infinite number of additional gluons, we still refer to this 
resummation as tree- level, since one resums leading combinations of coupling constant and field (gA). Although their 
momentum is integrated over, these gluons don't appear in loops but as additional gluons connecting the soft and 
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hard parts and as such are at the same level as other partons. The procedure is in fact ensuring that the non-locality 
is color gauge invariant, i.e. the replacement id^ — > iD^ — id^ + gA^. The importance of the directions of gauge links 
was realized in Ref. Even in the case of a light-like non- locality, the links are either past- or future-pointing, but in 
the absence of transverse separation this feature becomes irrelevant in the squared amplitude. In the case of transverse 
momentum dependent (TMD) correlators |5[ and correspondingly TMD soft functions [y, [7j, one also must account 
for a transverse non-locality, requiring more complicated Wilson lines. These Wilson lines have been extensively 
studied The TMDs lead to a rich phenomenology of azimuthal asymmetries 0, 0, Q EH ■ The distinction 

of past- and future-pointing gauge links provides a natural explanation of single spin asymmetries at the partonic 
level. Through the gauge links, time-reversal odd (T-odd) parts are incorporated in the TMD correlators and the soft 
functions in their parametrization within a field theoretical framework of Quantum Chromodynamics (QCD). Beyond 
tree-level, many complications arise [l6l - fl8j , of which at present certainly not all implications have been investigated. 
Although higher orders in QCD may invalidate any tree-level results, we will follow here the diagrammatic approach 
outlined in the above that provides us with the basic field theoretical picture which needs verification in an all-order 
QCD treatment. 

In section [HI we introduce some of the basics of TMDs needed in the rest of the paper. As said, a particularly 
interesting feature for the TMD soft functions entering the description of hard processes already at tree-level is the 
non-trivial nature of the Wilson lines connecting the non-local field operators in the matrix element. This structure of 
Wilson lines, arising from both collinearly and some of the transversely polarized gluon fields, now becomes dependent 
on the color flow in the hard process. In particular, in those situations that multiple color flow possibilities exist, 
this gives rise to an entanglement that can spoil factorization already at tree-level. How one gets the basic tree-level 
entangled result in a diagrammatic approach is outlined in section [TTTl 

The aim of this paper is to show that this entanglement simplifies for the 1-parton unintegrated (1PU) case, by 
which we refer to a situation in which only the transverse momentum in one of the hadrons is manifest. What remains 
is at tree-level a factorized expression with a correlator that still does have process dependence, which is contained 
in a non-trivial process-dependent gauge link. This is sometimes referred to as generalized factorization. The proof 
is given in section IIVI We will also show how one can proceed if transverse parton momenta in several hadrons are 
involved. In such cases one can consider weighted asymmetries expressed in terms of transverse moments of the TMD 
functions. These weighted asymmetries actually again involve only collinear functions, but these functions are given by 
matrix elements of higher twist operators, among them gluonic pole matrix elements or Efremov-Teryaev-Qiu-Sterman 
(ETQS) functions. In section [Vl we will use transverse moments to analyse single and double weighted asymmetries 
and illustrate this for Drell-Yan (DY) scattering and for a process with quark-quark scattering as underlying hard 
partonic process. The complications and non-universality always involves the gluonic pole matrix elements, which 
have been extensively studied [19|-|26j. Since such matrix elements vanish for fragmentation [27H29j . we will not have 
to worry about the transverse momentum in the final state. We will comment on this further in our conclusions. 



II. COLLINEAR AND TRANSVERSE MOMENTUM DEPENDENT (TMD) CORRELATORS 

We have split up this paper in a number of sections, in which we discuss in a diagrammatic expansion the inclusion 
of all gluon fields that contribute at leading order and tree-level. The starting point is a hard subprocess, for which 
we will consider as the most generic example a two to two process with a truncated amplitude ^{px,P2\ ki, fo), from 
which the wave functions of the partons (Dirac spinors u(p\) for quarks, or polarizations e(pi) for gluons), are omitted. 
Rather than through wave functions, the external partons are accounted for through quark or gluon correlation or 
spectral functions, which are built from matrix elements of the form (X\tp(^)\P) involving hadron states \P) rather 
than a free parton wave function (0\ijj(t;)\p). This immediately brings in the need to also consider multi-parton matrix 
elements with the same states, such as (X\A fJ '(ri) tp(^)\P). 

These matrix elements appear as squared contributions in the correlators (including Dirac space indices i and j), 

= eip< (p& 3 (0)MO\p), (i) 

pictorially represented in Fig. [lja). Usually, a summation over color indices is understood. This means that we will 
have $(p) = Tr c [$(p)], where $ij(p) is considered also a matrix in color space, made explicit $jj ;rs oc Vir(£) Vj S (0)- 
Including gluon fields one has quark-quark-gluon correlators like 

^ Aij (p, Pl ;P) = J d^dS e * e*^ (P|^(0) A"(r,) M0\ p ), (2) 
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illustrated in Fig.QJb), and similarly matrix elements with more partons. The color structure of the field combination 
ipr^ip s (0) in the quark-quark-gluon correlator now actually has a color octet structure, denoted (when appropriate) 
as $ 8 = Tr c [<£> T a ]T a . Using for = A^ a T a a matrix-valued field we have $^ = Tr c [$ 8 A^]. In some cases, 
it will be convenient to explicitly use the momentum space fields defined as ip{p) = J d A £, e zp '^ ^(0 an d A«(p) 
| <j 4 £ e lp '^ ^4 M (£), which in the case of free fields would have parton and anti-parton contributions multiplying on-shell 
factors (27r) S(p 2 — to 2 ) 6(±p°). For the correlators, we get 

{2n)H\p- P ') ZufaP) = (p$ 3 (p>)Mp)\p), (3) 

{2-Kf5\p-p')^ Aij {p, Pl -P) = J ±^(P\4i j {pi)A»{p 1 )Up-Pi)\P)- (4) 
The corresponding correlator describing fragmentation into hadrons is for quarks given by 
Ay(k;K h ) = E(^)I / e^ k < (0\M0)\Kh,X)(K h7 X\^(O\0) 

= (2^ J d ^ C ^ k< (°l^(0)4^(0|0), (5) 

pictorially represented by the blob in Fig.QJc). An averaging over color indices is implicit, thus A(fc) = -^-Tr c [A(fc)] 
with again A(fc) a diagonal matrix in color space. The second expression in the above involves hadronic creation and 
annihilation operators ajjO) = \Kh). In a momentum space representation for the operators, we have 

(2n) A 5 i (k - k') Aij (k; K h ) = ^(0|^(fc)|Jf h , X)(K h , X|^(fc')|0). (6) 

In fragmentation correlators, one no longer deals with plane wave hadronic states, but with out-states \Kh,X). 

We want to give an expression for the cross section of the (semi)-inclusive process H\(P\) + ^2(^2) — > h\{K\) + 
^2(^2) + • • • in a kinematic regime where PfPi ~ P%-P2 ~ K\-k\ ~ K^-ki are small (we will refer to this scale as 
the squared hadronic mass scale, M 2 ~ 1 GeV 2 ) as compared to the usual hard invariants in the full or the partonic 
process such as s ~ 2P1/P2, t\ « —2K\-P\, s k, 2pi-p2, t w — 2ky P i (we will refer to this scale as the squared hard 
scale, Q 2 3> M 2 ). We note in passing that in cases where heavy quarks are involved, those quark masses of course have 
to be included. In a hard process as described here, we aim for a description in which the squared partonic amplitude 
p#| 2 is convoluted with the correlators $(p, P), A(fc, Kh), etc. In order to illustrate the use of the correlators, assume 
an amplitude ^# cx Fi T 2 as illustrated in Fig. [3Ja). The simplest tree-level diagrammatic contributions to the cross 
section that can be written down is shown in Fig. Hfb) and is of the form 

d<7~Tr c [$(p 1 )rtA(fc 1 )r 1 ] Tr c [$(p2)r^A(fc2)r 2 ], (7) 

where Tr c [. . .] parts are traced over color. This expression still needs to be integrated over the parton momenta, 
which will be discussed below. In the case that the vertices T don't have any color structure, one can, because of 
the simple color singlet structure of $ and A in the quark-quark correlators, perform the color trace separately for i> 
and A, Tr c [$(p) T* A(fc) V] = Tr c [$(p)]^-Tr c [A(fc)] V V* (one summed and one averaged) and the cross section can 
be written in terms of the color-traced entities 

da ~ $(pi) *(pa) r x rj r 2 r* A(h) A(k 2 ), (8) 
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FIG. 1: The pictorial momentum space representation of quark-quark correlators for distribution functions (a), a quark-quark- 
gluon correlator (b) and a quark-quark correlator for fragmentation functions (c). 
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where the remaining contractions are Dirac space and Lorentz indices, which have been suppressed in both Eqs[7]and 
ISJ These expressions of course should be extended with all possible correlators containing quark and gluon fields, in 
which cases color traces become more complicated because the quark-quark part of the correlator can have a color 
octet structure ($s)- The restriction to hard kinematics limits the number of diagrammatic contributions, although 
even at leading order, there still are many gluon contributions as will be discussed in Section ITiil These will complicate 
the color tracing, an effect which is particularly important when dependence on transverse momenta is considered, 
which precisely is our goal in this paper. We will only work out leading contributions in an M/Q expansion, although 
the separation of various orders requires care, as we will argue in Section [TTTl The M 2 /Q 2 effects certainly cannot be 
calculated in our diagrammatic approach. At that level, there are many contributions that spoil already at tree-level 
the possibility to write down in a consistent way a cross section in the form of Eq. and certainly to perform the 
color traces as in Eq. [3] 

For parton momenta relevant in a hadron correlator (hadron momentum P) we make the Sudakov decomposition, 



p = x P + p T + (p-P -xM 2 ) n, 



(9) 



where the role of the (approximately) light-like vector n, satisfying P-n = 1 can come from any of the hard (external) 
momenta, e.g. n = K^jK^-P or n = k/k-P (provided k-P ~ K^-P ~ Q 2 )- The momentum fraction x = p-n = p n is 
^(1). For any contractions with vectors outside the correlator $(p, P) one has P ~ Q, p T ~ M and n ~ 1/Q. Note 
that if n is an exact light-like vector, one can construct two exact conjugate null-vectors, 



= P — A M 2 



and 



(10) 



satifying n + -n_ = 1 and ni_ = n_ = 0, that can be used to define light-cone components a = a-n^i (thus x = p-n 
p + ). The symmetric and antisymmetric 'transverse' projectors are defined as 



'? = - nfn] = g" u - P^n !/ > + M 2 n^n v » g" u - pt"n">, 



(ii) 

(12) 



Since transverse momentum dependence is a central issue in this paper, we have to worry about different n vectors. 
With Art = n' — n ~ 1/Q one has at 0{Q°) that Aa; w Ap 2 ^ w although the transverse momentum itself does 
change at order 0(1), Ap T = —Ax P. 
The integration over parton momenta, 



P 



dx d p T da 



d(p-n) d 2 p T d(p-P), 



(13) 



is insensitive to the particular n vector. In view of the relative importance of the components in this integration, one 
can, upon neglecting any M 2 /Q 2 contributions in the cross section, integrate within a soft correlator over p-P (i.e. 





(a) 



(b) 



FIG. 2: For the purpose of illustrating the structure of Wilson lines, we use the hard amplitude with one particular color flow 
for the quark lines as shown in (a). The squared amplitude needed for the cross section of the scattering process initiated by 
two hadrons with momenta Pi and P2 is shown in (b). 
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p ) to obtain the TMD correlator 

$(x,p T ;n) = J dp-P $(p;P) = | e'"* (P|^(0) ^»(0|P> 



(14) 



which we will still consider as the unintegrated correlator. On the left-hand side the dependence on hadron momentum 
P has been suppressed. In the TMD correlator the non-locality is restricted to the light-front (LF: £-n = £ + = 0) 
and the correlator depends oni= p-n and p T . The light-cone correlators are the collinear correlators containing the 
parton distribution functions depending only on the light-cone momentum fraction x, obtained upon integration over 
both p-P and p T , 



$(x; n) 



dp-Pd 2 p T $(p;P) 



(27T) 



(PlV(o)V(Oli') 



(15) 



LC 



where the subscript LC refers to light-cone, implying £-n = £ T = 0. This integration is generally allowed in hard 
processes up to M 2 /Q 2 contributions and also up to contributions coming from the tails, e.g. logarithmically divergent 
contributions proportional to cc s (Pt)/Pt tails (30j and relevant when looking at evolution |31| . Such contributions, 
however, require next-to-leading order (NLO) QCD, which goes beyond the tree-level resummations that we discuss 
in this paper. In diagrammatic language they for instance involve ladder graphs describing emission of gluons into 
the final state, relevant for the evolution of the correlators. The collinear correlators are relevant in hard processes 
in which only hard scales (large invariants ~ Q 2 or ratios thereof, angles, rapidities) are measured. If one considers 
hadronic scale observables (correlations or transverse momenta in jets, slightly off-collinear configurations) one will 
need the TMD correlators for a full treatment. 

The correlators encompass the information on the soft parts. They depend on the hadron and quark momenta 
P and p (and in general also spin vectors). Depending on the Lorentz and Dirac structure of the matrix elements 
involved one can look for the pieces in the correlator that show up as the most dominant matrix elements among 
the contributions in the hard process. Including also gluon fields, the Fourier transform of matrix elements with a 
maximal number of contractions with n, 



<^(0)#(0> and ( G""(0)G^(0 ), 



(16) 



(the latter with transverse indices a and 0) are the dominant combinations that appear in the correlators. They 
are the dominant ones because the contractions with n lower the canonical dimension of the operator combination, 
minimizing the power of M that after contractions of open indices inevitably is the scale of the hadronic matrix 
elements. The two matrix elements above have canonical dimension two. The corresponding local matrix elements, 
ip(0) i/ttp(0) and G na (Q) G n/3 (0) for quarks and gluons, respectively, are color gauge-invariant (twist 2) operators, the 
non-local combinations in Eq. [Tulare not gauge invariant. Expanded into local operators, the expansion would involve 
operator combinations with derivatives such as ip(0) ^ d n ... d n ip(0). Color gauge invariance in the correlators requires 
in the local matrix elements covariant derivatives or in the non-local matrix elements the presence of a gauge link 
connecting the two fields. For the light-cone correlators the gauge link corresponds to the inclusion of an arbitrary 
number of 'leading' gluon fields A n (rf) in the field combinations in Eq. \W\ which are resummed into a gauge link 

m m = Tr c [W™ j mm] , given by 



= 9 exp 



i / drj-P n-A(r]) 



(17) 



Including this gauge link, the non-local operator combinations 



and 



LC 



K,o] 



LC 



(18) 



can be expanded into twist two operators "0(0) i/iD n ■ ■ ■ D n ip{0) and G na (0) D n . . . D n G n/3 (0) for quarks and gluons, 
respectively (number of D n, s is the spin of these operators). Also TMD correlators require a gauge link, but the 
separation of the two fields is no longer a simple light-like one and they involve derivatives with transverse indices. It 
is important to realize that in principle any gauge link with an arbitrary path gives a gauge-invariant combination. 
What is the appropriate link contributing at leading order (in M/Q) in a given hard scattering process, however, is 
calculable (see next section). 

The correlators in this section have been rewritten into matrix elements of non-local products of fields. They involve 
both quark and gluon fields as well as hadronic states or hadronic creation and annihilation operators. For them there 
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FIG. 3: Inclusion of collinear gluons from $a...a(p — pi . . . — Pn,Pi, ■ ■ ■ ,Pn) coupling to an outgoing (colored) quark line with 
momentum k. 

is no systematic perturbative expansion in terms of the strong coupling constant. The kinematic separation of soft and 
hard, however, allows the integration over p~ = p-P, leaving a product of operators at the light-front, that is at equal 
light-cone time £ + = £-n = 0. For such a product the time-ordering is automatic, which means that the p~-integrated 
parton correlators, thus, can be considered as a cut anti-parton-hadron scattering amplitude, i.e. a Green function, 
untruncated in the parton legs [32| . This is the case for both collinear and TMD correlators (33[. This identification 
has been very important in deep inelastic processes [34j . allowing the use of analyticity and unitarity properties of 
field theories, at least under the assumption that these properties apply to QCD. We will use it later in this paper for 
fragmentation correlators. 



III. COLOR GAUGE INVARIANCE 



It is worthwhile to repeat the steps that lead to color gauge invariance by including collinear gluon fields, which 
will also be the first step to obtain the Wilson lines in TMD correlators. We include in the diagrammatic approach 
matrix elements with gluon fields for which we make a Sudakov expansion, 



A"(r]) = A n (n) P" + A^rj) + (A p (n) - A n (rj) M 2 ) n» 



(19) 



A similar expansion can be written down for A fJ, (p). In order to use Ward identities it will be convenient to look at 
the (collinear) gluon field component along parton momentum p^, hence we write 



d 4 rj e 



ip-ii 



A n (r]) 
p-n 



(p-n) A% (n) - K A n (rj) {p-n) A p (n) - (p-P) A n (n) 



p-n 



p-n 



In the correlator the momentum p^ 



id^^q), so 



A^ip) = — [ a 

p-n J 



4 n e tp - v 



1 

p-n 



A n (n)p tl + id n (n) A%(n) - id^(n) A n (n) + (id n (n)A p {n) - id F '(rj)A n (ri)) 
A n {p) p^ + i G^(p) + i G nP ( P ) n» 



(20) 

(21) 
(22) 



Although the latter appears to be only true for the Abelian case, we will find the same result in the non-abelian 
case, but to complete that proof, we first need to incorporate the collinear gluons into the matrix elements. Using 
the expansion in Eq. [5U rather than the one in Eq. [TO] streamlines the inclusion of collinear gluons circumventing the 
explicit treatment of transverse momentum dependent parts (as done in Ref. @). The results are of course identical. 



A. Collinear gluons 



Since matrix elements involving ip(0) A n (r}i) . . . A n (n N ) ?/>(£) and G" Q (0) A n (r)i) . . . A n (n N ) G n0 (£) for quarks and 
gluons are as leading as the matrix elements without collinear gluons, the contributions of correlators with gluon fields 
of which the polarization is along its momentum p^, i.e. the first term in Eq. 1221 need to be included in the leading 
expression for the cross section. Using Ward identities, only the contributions of gluons coupling to external parton 
lines survive as discussed in detail in Ref. [Tlj . The resummation is best illustrated by looking at the example of 
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gluons attached to a correlator attached to an outgoing quark line with momentum k (see Fig. [3]). The results 
of these parts within the amplitude are 

A = ^(fc)r^(p), 

where the ip(p) and ^>(fc) are fields belonging to the correlators of initial (momentum p) and final state quark (mo- 
mentum k) respectively. The one-gluon contribution to the link is 

J (2tt) 4 pi-k {-2k-pi+ie) 

The numerator becomes — ^i) = = {$,$1} = 2k-p\. The added term is zero since ip(k)$ w 0. Thus one has 
(note that the sign of k ■ pi is positive) , 

Ai = m I , f iPl l r^(p- Pl ) = m [/j. fc](1) TMp), (23) 
V (27r) 4 (-k-pi+ie) lI 1 

which is the one-gluon contribution to a path-ordered exponential (see Appendix [XJ. The two-gluon term becomes 

J (27r) 4 (27r) 4 pi ■ k P2 ' k ft -pi f-pi-fc 

d 4 P i d±p 2 - A k { Pl ) A k ( P2 ) hAtljK) Mt-h-fc) Tl( _ v 
(2tt) 4 (2tt) 4 V{ ' Pl ■ k p 2 -k {-2k ■ Pl + te) (-2k ■ Pl - 2k ■ p 2 + ie) WKP Pl P2> 

- f dV d 4 P2 A k ( Pl ) A k (p 2 ) 

= m J W (2^ (-k-p 1+l e) (-k.p 1 -k. P2+ *e) T ^ p - p ^ P2) 

= ^(k) U l + ]{2) Tijj (p). (24) 

This is the two-gluon contribution to the path-ordered exponential. This term illustrates the recursive procedure that 
gives the all-gluon results to the full path-ordered exponential or Wilson line (see Appendix |A|) . 

00 

Y,A N =mU [ * ] riP(p), (25) 

N=0 

where 

U [ + ] ijj(p) = J d*f exp (ip ■ £) ^exp ^-ig ^ d(r) ■ P) A™ (77)^ (26) 

In] 

is the Fourier transform of the field including the Wilson line of which the indices of U + indicate its direction, namely 
a Wilson line running from light-cone +00 to £ along the light-like direction n. We will instead of the notation with 
the contraction, mostly use the notation 

~$(k) U [ + ] Tj; {p) = ^(k) U l + ] [p] Ttp(p), (27) 

The second expression gives the Wilson line a label [p], indicating that all fields in it belong to the correlator $(p), to 
which also tp(p) belongs. The problem with this path-ordered exponential is that it is a unitary matrix in color space 
that is stuck in the respective traces, although it is in fact only the (symmetric) color charge operators T ai . . . T 0,3 of 

the terms JJ+^ in the expansion of the U+ that are stuck there. 

For the gluon insertions coming from a particular correlator and coupling to an incoming fermion line one finds a 
Wilson line connecting to light-cone —00, which is a consequence of the sign k ■ pi being negative in that situation. 
Including all multi-gluon interactions originating from Q(pi) in Fig. [2f b) , we get the diagrammatic result 



da ~ , & c [*(pi)nc4 fcllt b'i]A(*;i)c4 fcl] |pi]ri] 



x Tr c [U [ ^[ Pl ] <S>(p 2 ) [px] T* U [ ^[ Pl ] A(k 2 ) U [ ^ ] \p x ] T 2 ] , (28) 

with the (color charge of the) Wilson line stuck in the color traces at the 'positions' corresponding to the external 
parton lines. Note that in Eq. [3E] the Wilson lines in sub-expressions like U^ 1 ^ [pi]A(ki)U+^ [p{\ are part of the 
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(a) (b) 

FIG. 4: (a) An example of a diagram with two gluons attaching to the same outgoing line with momentum k^. (b) In more detail 
these couplings are discussed in the appendices and they result into one gauge connection I/!? 2 ' [pi, P2, fci], which combines all 
collinear gluons coming from < I > (pi), ^(jpi) and A(fci). Together with transverse pieces, and combining it with gauge connections 
at other legs one will get a full color gauge invariant result, discussed in Section Hill 

Fourier transform in the correlator In coordinate space, taking £1 to be the coordinate conjugate to pi, the 

link and conjugate link actually run between different points, in this case [/+ in U^t 1 ^ [pi] ■ ■ ■ ip(pi) corresponds actually 
to f/r*^, ti, while the conjugate link U+ in tp(p\) . . . [p{\ corresponds to uIq^ +00 y Wilson lines between 0i T and 
£i T are still lacking, so Eq. [28] is certainly not color gauge-invariant. 

The next step in our treatment is the inclusion of gluon interactions coming from different correlators, say $(pi) 
and <I>(p 2 ), coupling to the same quark line with momentum k 2 , such as shown e.g. in Fig. SJa). They give rise to 
intertwined Wilson lines. Two examples of this have been given in Appcndix[B] for the case of two gluons coming from 
different correlators and for the case of three gluons, two coming from the same correlator and one from a different 
correlator. These examples illustrate the recursive procedure. The result for all insertions to a particular leg is a color 
symmetric combination of the insertions from all correlators, which since all gluon polarizations are identical is just 

U^lp^h] = SiU^lp^U^lpip^lh}}, (29) 

in which the ordering of the three connections on the right-hand side is irrelevant (fully symmetrized). This is 
illustrated in Fig. 0] Including all multi-gluon interactions from <J>(pi), A(A;i), $(p 2 ) and A(fc 2 ) onto all legs, Eq. |2"51 
generalizes to 

da ~ Tr c [C/I pilt [p 2 , h,k 2 ] <t>( Pl ) U M \p 2 , k u k 2 ] rj U^[ Pl ,p 2 , k 2 }A(kt) U^ l] \p u p 2 , k 2 ] I\] 

x Tr c [U^^MM HP2) U M \jn,kiM n U [ ^[ Pl ,p 2 , fcj A(fc 2 ) U^ ] [ Pl ,p 2 , fcj T 2 ] , (30) 

illustrated in Fig. EJa). We note that the Wilson lines in Eq. l30l have different light-like directions, which originate 
from the fact that we in the decomposition of gluon fields in Eq. l22l simplv made the most convenient choice depending 
on the particular correlator. For two different Sudakov decompositions of the gluon field, one finds that 

A n (p) _ A n '(p) = iG nn '(p) 
p-n p-n' {p-n)(p-n') 

The field G nn (p), however, appearing in a correlator will not contribute at leading order, but at subleading 

order (1/Q). This allows one to replace all the light-like direction dependence in Eq. [3D] by a generic null-vector n. 
Even after choosing one light-like direction, the result in Eq. [30] is not yet gauge-invariant, since transverse gauge 
connections are still missing. 

B. Collinear correlators 

The missing transverse pieces don't matter when one takes a collinear approach, implying integration over transverse 
momenta p T besides the integration over p ■ P. In that case there is no transverse separation of the fields. After 
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integration over transverse momenta, one has a color gauge-invariant result, 

a ~ Tr c [[/ W V, k u k 2 ] *(xi) U l l l] [p 2 , h, k 2 ] rj E/j nIt [pi,P2, felAfa) C/|" ] \puP2, k 2 ] T x ] 

x Tr c [C/I nI+ [pi, hi, k 2 ] $(x 2 ) f/ W [pi, fc x , k 2 ] T* E/f ]t [pi.P2> fci] A(z 2 ) U [ + ] \pi,p 2 , k{\ T 2 ] . 



(32) 



One can combine the Wilson lines to and from light-cone ±00, all made up of A n fields, into finite Wilson lines, e.g. 

= t r + l '^[pi] f+^Pi] since after the integration over p\ T they not only both run along n, but they coincide 
since one also has T = £ T . Furthermore, it is irrelevant if one composes from Wilson lines running via plus or 
via minus infinity, and also the direction n is in fact irrelevant. It is just the direction of the straight line connecting 
and £. We recall that the argument pi or 11, given to the Wilson lines, is simply needed to indicate that the fields 
in that Wilson line belong to the correlator which is the Fourier transform of the matrix element (?/>(0)V>(£)). 

Thus in coordinate space one just has the Wilson line in Eq.[17l which connects the points and £ in composed 
of WrE?L = f7rnL_i f7 r „ „. As far as relevant for $(xi), the Wilson lines in the first trace form a gauge link, those in 



[0,«] w [0,oo] "[«>,£]• 

the second trace form a closed loop, which in the collinear situation (when 3 
space. One is left with 



a ~ Tr c [U [ ^[k 1 \^(x 1 )U i y[k 1 }rlU i y[ Pl }A(z 1 )U l : il [p 1 }] 



Tr c [[/ Wt [fc 2 ]$(a; 2 )(7l 



[k 2 ] T* 2 U. 



„]t 



P2 ]A(z 2 )i7l nI [p 2 ]r 2 ] 



£ T ) becomes a unit operator in color 



(33) 



The way of turning the gauge connections into gauge links at the collinear stage is actually just applying gauge 
transformations U,„ (with a fixed point a) to all fields. This can actually directly be applied to Eq.[321 One obtains 



r £, (with a fixed point a) to all fields. This can actually directly be applied to Eq. 

a~$w\ Xl )<f>W( X2 ) r x rjr 2 r* ti w \ Zl )ti w \z 2 ), 



(34) 



where 



$W(x) = T* e [U£ ] \p]*(x)U^\p]] 



Tr c [VF w [p]$(a;)] 



d£,-P 
2tt 



d£-P 
2tt 

Jp< 1 



e^< (p\mu^ ±oc] uiii s m\p) 
p\mw^(o\p) 



LC 



(35) 



LC 



and 



N, 



Tr c [4" 1 [fc]A(z)4" lt [fc]] = ±-Ti c [A(z)wW[k]] 



dtj-K h 
2tt 



N, 



Tr c (0|^ l oo0] ^(0)aj l a^(Of/£ l ±oo] |0) 



(36) 



LC 



are the color gauge-invariant collinear correlators, including unique gauge links W along the light-like separation. The 
gauge link being unique, it is usually omitted. These (color gauge-invariant) correlators can be expanded in terms of 
the standard parton distribution functions and fragmentation functions, respectively. 

For the correlator $(x) it is also possible to circumvent manipulating with Wilson lines by taking moments in 
x\ = p\-n, x 2 — p 2 -n, 1/zi — k\-n and l/z 2 = k 2 -n. Using 



7 r) n W m 



Kin wn &> 



(37) 



one easily verifies the factorized expression in Eq. 1341 In the process of taking moments one then encounters x N (x) 
involving matrix elements with covariant derivatives D n , 



x N & w \x) 



d£_-P 
2tt 



•>« {p\mw™ G (iD n ) N M)\p) 



(38) 



LC 



We have shown everything for one correlator, but one has similar expressions for the other correlators. Since the 
gauge link is unique for the collinear correlators, we will just write $(a;i), <f>(x 2 ), A(zi) and A(z 2 ). 
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0(6) : ^(0 2 ) 




FIG. 5: (a) The gauge connections from all collinear gluons from the various soft correlators produce (entangled) gauge 
connections with color charges located on the external legs of the hard part, (b) We have indicated for the correlators and 
gauge connections also the actual space-time points they are bridging limiting ourselves for simplicity to the coordinates 
conjugate to pi (points 0i and £i; see also the discussion following Eq. l28|) and p2 (points O2 and £2), leaving out the space-time 
structure for the fragmentation correlators for which we would have to include also the coordinates conjugate to fci and ki. 



C. Including transverse gauge connections 



To see how transverse gauge connections arise, we jump back to Eq. 1301 and note that not only the (. . . "0(0) ^KO • ■ ■) 
correlator has acquired gauge connections along the n-direction, but also other matrix elements involving the other 
gluon components in Eq. 1211 e.g. the operator combination 

(■■■i^W^A^U^ -id^ri)U^ v] A-{ n )U^ M ..), 

for which we can use the non-abelian relation (directly based on Eq. [37| , 



(39) 



Differentiation with respect to 9" = djdrf gives 

(. ..i^)^ A?® -i^U^ A n (77) [/W 



showing that the transition of Eq. [51] to Eq. [22] also works in the non-abelian case, if the appropriate Wilson lines are 
included. Similarly, one has 



p 

rfirt AP w <L, - ^V =00 = £ dc p < C] (g- p (o + [d p , A«( C )]) <L] 



(40) 



C, n =n n Xt=vt 



The subtractions at light-cone infinity are important to produce the missing transverse Wilson lines. It is shown 
in detail in Ref. [9J how these subtractions A"(i])\ rl p =00 from the 'higher twist' matrix elements with A T (and A p ) 
fields provide the missing transverse gauge connections U^ T ^ = U^ t U^ OQt ^ 1, which we will indicate as U T [p] 
with the argument p again just indicating that the link connection involves the (transverse) endpoints T and £ T 
conjugate to parton momentum p and fields belonging to &(p). Incorporating all transverse gauge connections, we 
obtain (symmetrized) pieces 



tr>i,p 2 , h] = s{u T [ Pl ] u T [ P2 ] u T [h}}, 



(41) 



with irrelevant ordering among themselves. In the ordering of these pieces and the combination with collinear pieces 
of the gauge link, we will neglect contributions from commutators involving A fJ '(pi) and A p (P2) since (pi —P2) 2 ~ Q 2 ■ 
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FIG. 6: The gauge connections present in the expression for the unintegrated cross section in Eq. 321 (a) U\ U+ => Wj™-' = 
<L, t/[o T , 5T] (b) C/l f/- => ^ Wt = ^,o Tl ^L,o]. (c) = W W W 



[» li- 



lt is clear that in kinematic regions where all parton momenta are small this cannot be used, complicating the full 
QCD analysis. Combined with the connections along the light-like direction, we now have for our tree-level analysis 
a set of entangled Wilson lines that bridge the non- locality of parton fields in the various correlators (see Fig- E](b)). 
The resulting expression for the cross section is again of the form in Equation Eq . 1301 but including transverse pieces, 
denoted U + [p, . . .] — U T [p, . . .] W n ' [p, . . .], etc. One gets for the cross section an expression in terms of unintegrated 
correlators, which are only integrated over p-P, 

da - Tr c [UL[p 2l fa,k 2 ] $(x 1: p lT ) U-{p 2 , fa,k 2 ] T* U\\pi,p 2 , k 2 ]A(zi, k lT ) U+\pi,p 2 , k 2 ] rj 

x Tr c [ul [p! , fci, k 2 ] #(x 2 , P2t ) V- [pi ,ki,k 2 ]T* 2 u\ [pi, p 2 , h] A(z 2 , k 2T ) U+ [ Pl , p 2 , k ± ] T 2 ] . (42) 

This resulting expression is now color gauge-invariant. The Wilson lines can be taken along a generic n-direction, 
which even could be chosen different for each of the gauge connections in Eq. |3H but its color structure is fully 
entangled and it does not allow for a factorized expression with universal correlators that have their own gauge links. 
Viewing it as a factorized expression it contains hard amplitudes, soft correlators and gauge connections, where the 
gauge connections take care of a 'color resetting' which feels all hadrons that are involved. 



IV. DISENTANGLING THE COLOR FLOW DEPENDENCE 



Starting with the expression of Eq. for the cross section, we first look at the case that the hard process doesn't 
affect the color flow as in Fig. [^a). This means that the vertices Ti and F2 are color independent. In our expression 
we have floating around in the full expression pieces, like U T [p] U+\p] . . . $(x,p T ) . . . E/+ [p] C/ T ^[p] ■ • ■, which only 
when they are combined in a single trace would yield gauge links appropriate for TMDs, 



<&M (x,p T ) = Tr c [U T [p] U l ; l] [p] p T ) U [ f[p] U T ^ ]p]] 

<^0^ (Pim u^ T] f/M s] mP) 



(2tt) 3 



LF 

Tr c [^i" 1 [p]$( 2 ;,p T ) 



(43) 
(44) 



LF 



or 



AW(z,fc T ) = ^Tr c [U T [k]U^\k]A(z,k T )U^[k]U^[k]} = ^-Tv c [A(z,k T )W^[k] 



d^-K h d 2 j T 
(2tt) 3 



lk< W ^(0\U^ TfiT] Ulll t0] ma{a h m Ug ±oo] U[ 6t>00t] |0) 



(45) 



LF 



with the coordinate space structure of W±L ^ shown in Figs^a) and (b). Even if one is in the lucky situation that 
one can combine the relevant Wilson lines in a single trace, there still would be combinations left in other color traces 
which form Wilson loops = wj™ VfI™^ as illustrated in Fig. HJc). 



A. TMDs in 'elementary' processes 



Considering, as a reference, first the 'elementary' Drell-Yan (DY) process, which like other electroweak processes 
such as leptoproduction or electron-positron annihilation, is relatively simple, because the color flow through initial 
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FIG. 7: (a) The color-flow in the Drell-Yan process and (b) the gauge connections including their coordinate space structure. 



and final states consists of a single color loop. In that case, taking the DY process as an example (Fig. [7]), the cross 
section is given by 



da DY ~ Tr c [C/ Wt [p 2 ] <f(x uPlT ) U m [p 2 ] V* U ml [ Pl ] <S>(x 2 ,p 2T ) U m [ Pl ] V 

= Tv c [w l ^[ P2 ] Hx uPlT ) r*$(x 2 , P2T ) w H [ Pl ] r] . 



"]tr 



(46) 



The gauge connections contain collinear and transverse fields and one has to be careful in disentangling. We note in 
particular that in general 



da DY + Tr c [lTi"V]*(zi,:PiT)] Tr c [$(z 2) p 2r ) IT Nt [p 2 ]] ^ TT 



(47) 



However, if one of the hadrons can be treated collinear, i.e. integrating over one of the transverse momenta, say p 2T , 
we can disentangle the traces and obtain 

d(TDY 



$ [ \xi,p lT )<S>{x 2 )Y, DY , 



(48) 



where $' _ '(pi) = Tr c [pi] $(pi)] . This illustrates the TMD factorization at tree-level and the relevant gauge 
links to be used, which was considered in detail in Ref. (35[. Beyond tree- level one has the factorization formalism of 
Collins, Soper and Sterman (36j, which has proven to be very successful for many applications, but does not catch all 
subtleties at small transverse momenta [l8| . We will return to the problems with finding a full tree-level factorization 
in two TMDs in Section [V] 



B. TMDs in 1-parton unintegrated processes 



In this subsection, we will turn to the situation of a more complicated color flow and we will show that TMDs 
with non-trivial gauge links appear in 1-parton unintegrated (1PU) contributions to the cross section, by which we 
mean contributions in which one has integrated over the parton transverse momenta in all but one of the hadrons. 
At this point we do not worry about experimental feasibility of measuring such a cross section. First assume as an 
intermediate step that only the transverse momentum of correlators $(pi) and $(^2) are left unintegrated, i.e. we 
look at jet production or non-hadronic final states. In that case one integrates over k\ T and k 2T and the collinear 
correlators A(zi) and A(z 2 ) just involve unique collinear gauge links, while Wb[fci] and VFrj^] become unity (in 
color space). The result is an expression of the form 

da! - Tr c [t/![p 2 ]$(a; 1 ,p lr );7_[p 2 ]rj;7;[p 1 ,p 2 ]A[ M/ ](z 1 )[/ + [p 1 ,p 2 ]r 1 ] 

x Tr c [£/i[pi] $(a; 2 ,p 2r ) E7_[pi] r 2 u\ [ P i,p 2 ] A^(z 2 ) U+[p u p 2 ] T 2 ) , (49) 
which is still entangled. This was the example discussed in Ref. [371 ] . 



13 



FIG. 8: The squared amplitude for the hard interaction between 
two quarks coming from a one-gluon-exchange diagram. We 
consider for simplicity two quarks with different flavors. 



Integrating over p 2T , the result reduces to 

da l ~ Tr c [$( Xl , PlT )rj ul[ Pl ]^ w \ Zl ) U+[ Pl ] T,] Tt c [u1_[ Pi ] & w \x 2 ) E/_[pi] ^ U^} A^ w \z 2 ) U+[ Pl ] T 2 ] 

~ TT c [<f(x uPlT )Tl AW( Zl ) W+[ Pl ) T 1 ] Tv c [wl[ Pl ] ^(x 2 )T* 2 A^ w \z 2 ) W+[ Pl ]T 2 ] 

= Tr c [W [ + ] [px] $(x uPlT ) rj A(zi) r x ] Tr c [W^ 1 [p x ] $(z 2 ) V* A(z 2 ) T 2 ] , (50) 

and gives 

da l ~ Tr^^^^i^^^inA^Orx] 1 ^Tr c [^ 1 [p 1 ]]Tr c [$( a;2 )r;A(z 2 )r 2 ] (51) 

= ^^(x^)^) TiT^T* A(zx)A(z 2 ), (52) 

S v ' 

Si 

where Wilson lines combined with the correlators can be combined into TMD correlators with gauge links, such as 
( P i) = Tr c [$(pi) [pi]] . In the step from l50l to ISTl we have split the second trace into two parts containing the 

trace of (appropriately averaged) and the trace of the (n-independent) collinear part. In the last step the trace 
of the Wilson loop (non-local in the coordinate £ conjugate to the parton momentum pi) has been absorbed into a 
more complicated gauge link for $[+( n )] (x x , pi T ) = Tr c [wj^ [ P i]] Tr c [$(pi) W+™' [ P i]] . Note that, in spite of what 

the notation might suggest, is not a multiplicative factor times because the fields in both factors belong 

to the same matrix element (indicated through argument [ P i]). These complex gauge links are the process-dependent 
gauge links discussed in Refs [Iol - [l2l \3§L |39| . 

As a further illustration, we include the color flow possibilities in the case of a hard process where a gluon is 
exchanged in the hard part (see Fig. [5]) instead of a colorless boson as in Fig. [5] One gets different color flow patterns 
that contribute giving rise to two different color contractions 

da ~ Tr c [£/ W V, k u k 2 ] <S>( Pl ) U M \p 2 , h,k 2 ] T a U l ^[ Pl , P2 , k 2 ] A(h) U [ + ] [ Pl , P2 , k 2 ] T 2 T b ] 

Tr c [C/ Wt [ Pl , kxM *(pa) U M [ Pl ,k u k 2 ] T* 2 T a U [ ^\px,p 2 , h] A(k 2 ) \p x ,p 2 , k{\ T 1 T b ] (53) 



X 



Nt+1 



Tr c [C/ Wt [ P2 , kxM *(Pl) U H [ P2 , kxM rj U l ^[ Pl , P2 , k 2 ] A(Ai) U [ ^\p u p 2 , k 2 ] T 2 ] 



x Tr c [C/ Wt [ Pl , hi, k 2 ] <5>(p 2 ) U^lpx, ki,k 2 ] T* U l ; i]f [ Pl:P2: k x ] A(fca) U [ ^[ Pl , P2 ,k{\ T x ] 
i Tr c [t/ Wt [p 2 , h,k 2 ] $(pi) C/ W [p 2 , k u k 2 ] VI U^\p uP2 , k 2 ] A(fci) U£ ] [pi,p 2 , fc 2 ] 



x Ut ]i \pi, fci, fca] $(p 2 ) (7 W [pi,fci, fc 2 ] T* (7|" ]t [pi, P2 , fci] A(fc 2 ) [/|" ] [p 1)P2) fcj] rj (54) 

= jf0i d<J i ~ 7ff=T da 2, (55) 

normalized by the usual color factor, such that the coefficients in front of these terms add to one. The first part is 
the one considered above. For the second part, one gets upon integration over k\ T and k 2T 

da 2 ~ Tr c [C/I" ]t [p 2 ] $(z 1)PlT ) U [ l l] [ P2 ] rj C/|" ]t bi,p 2 ] Afa) ^[pi,^] T 2 

x U™]pi] <S>(x 2 , P2t ) U H [ Pl ] r* ^[pi.pa] A(z 2 ) C/M^.Pa] I\] , (56) 

leaving an entangled situation. If only one transverse momentum remains unintegrated, one again can combine the 
result into a single complicated gauge link, 

da 2 ~ Tr c [$(xi,pi T ) rj C/i" ]t bi] A(*i) t^" 1 [pi] T 2 [/I^V] $(x 2 ) C/ W [p x ] T* t/j^bi] A(z 2 ) [pi] r x ] 

= Tr c [$(^ 1 ,p lT )^}" 1 [p 1 ]rj A(z 1 )r 2 $(x 2 )VF [ i ] " 1 [p 1 ]r;A(z 2 )r 1 ] (57) 
= & +a }(x 1 , PlT )$(x 2 )r 1 r* 1 r 2 r*A(z 1 )A{z 2 ), (58) 

s 2 
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where = Tr c [$(pi) wjT-'fpi] W^[pi]] . Combining both contributions one gets the result [Tlj . 

da ~ (^$ [+(D)I (*i,Pit) *(a; 2 ) ±1 - ^^(^uPir) *(a*) Sa) A(^) A(z 2 ). (59) 

In this particular case Si and £2 are the same, but in general one has more diagrams contributing to the hard process, 
each contribution gets split up into its own color patterns. The results for such 1PU processes have been tabulated 
in Ref. [n|. 



V. ANALYSIS IN TERMS OF TRANSVERSE MOMENTS 



In this section, we will start with the unintegrated result in Eq. 22] and take transverse moments, which will give 
collinear results, which include correlators with covariant derivatives iD" and A" or in gauge-invariant form integrals 
over G na . The transverse moments are obtained by looking at weighted cross sections of the form 

(p? 1 ...p?M = / d2 PT Pr 1 ■ • -Pr N -^~ t - (60) 

We will apply this to the cross section cr, built from the two pieces u\ and 02, discussed in the previous paragraph. 
The simplest, lowest, transverse moment is the integrated cross section, 

{a) = J d2pT = Hx2) ^ ^ -nTT ±2 \ A(zi) A(Z2) ' (61) 
which has between brackets the standard partonic cross section E. 



A. Single weighted asymmetries 

To see what happens with the first transverse moment, we need to investigate at the level of correlators what happens 
when one applies id^ to the matrix elements. To clarify the role of the gauge link in a correlator, in particular the 
contribution at infinity, we consider 

«S [■ • ■ <L] UprM • ■ • m = [■■■ <L] U[ 0t , (t] l D« . . . m] • (62) 

We can evaluate explicitly 

iu; u™ & . . . m = iD«(o . . . ^(o - f S d v .p u™ v] g™ {r,)u^ & . . . m, (63) 

J oo 

noting that the second term depends through the integration limits on the structure of the Wilson line, in particular 
on the transverse piece being at plus or minus infinity. We could simply use the second term as a definition of 
the A" field but such a definition would for instance not have the correct time-reversal property. We will use a 
time-reversal-definite expression for A a (^), 

/oo 
d V -P e (£ P - rf) U™ G na ( V )ul; ] ti] , (64) 
-oo 

which is gauge-equivalent to the expressions in Eqs 131] and [HE The form in Eq. [M] implicitly implies A" = and 
antisymmetric boundary conditions at £ p = ±00 for A"(£) and A p (£), 

u Us + 4too] A a (-oo, D^L,a = 0. ( 65 ) 

while 

/oo 
d^PG na (0 = 2nG na (^,C), (66) 
-OO 
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or combining these two equations, 



(67) 



The integral expression for TrG na represents actually the Fourier transform of the field strength tensor at zero mo- 
mentum in the n-component (which is indicated by including a tilde over the G and omitting £ p from the argument 
list) . Eq. becomes 

(68) 



To implement this in correlators, we consider besides the TMD correlator <$>(x,p T ) in Eq. 021 the collinear matrix 
elements <&'\{x, X\), which are the integrated quark-quark-gluon correlators starting with Eq. [2J Including Wilson 
lines, this collinear correlator is of the form 



d^-P drj-P 



2tt 2tt 



or integrated over Xi , 

$2(a;) = / ' dxx*'X(x,xi) 



,i (p-pi)-Z pipi-v 



d£-P 



LC 



— (p\mu^ A a (om\p) 



LC 



Similarly we define $>%(x, x\) and &q{x,x\) using iD a (rj) and G na (r]), respectively. The latter is given by 
d£-P df]-P 



$g(a;,a;i) 



2tt 2tt 

The relation in Eq. [M] implies for these correlators 

$1(X)=PV fd Xl —l—QZfax!). 
J X-Xx 



(69) 



(70) 



(71) 



(72) 



The operator G™(£ T ,£™) in Eq . appears in the correlator with a gluon field strength, but at zero n-componcnt 
of the momentum, i.e. at x\ — 0. The other (among them transverse) components of the gluon momentum are also 
integrated over in Eq. 1711 This particular correlator, referred to as a gluonic pole matrix element, 



d£-P 
2tt 



(73) 



LC 



will play an important role for transverse moments. Note that we have not included any argument for the field 
strenght in G na ip{x) , because all arguments become zero after integration over p\ T , p\-P (implying ^ T =r] T and 

LC 

£ > n =rj n ) followed by the integration over p T and p-P (putting £ on the light-cone, i.e. £ T =0 T and £ n =0). 

Using the expressions in Eqs and I63[ we consider the following unintegrated TMD functions starting with 
$W(x,p T ), 



$ a 9 l±] (x, PT ) = p%$M(x,p T ) 



(p\m u$ ±oo] u^ T] iD% uin^ m\p) 



LF 



Using the result in Eq. [68] one gets 



[±] (x, p T ) = ^ l±] (x, p T ) ± n^ [±] (x, p T ) , 



(74) 
(75) 



where 



^% [± \x,p T ) = ±J e<* W(0)^l tH ^ ] i4;(±oo,{ T )Er [ ^^(0| J P> 



d£-Pd 2 £ T 



LF 



:[±] 



(x,p T ) = 



d %p T ^ (p\mw [ ^d°m\p) 

d(,-Pd 2 £ 7 



LF 



LF 



JP< 



(2tt) 3 



{p\m w^, s] - ^(oip> 



(76) 



(77) 



LF 
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or after integration over p T , 

$^ [±1 (x) ee f d 2 p T p«$M(x,p T ) = <f>%(x)- I d Xl —-L-—*°( x , Xl ) (78) 



x — x\ =F it 

= $2 (z)±tt (79) 

S -~v ' 

*g(x) 

We note that the sole dependence on the direction of the gauge link is in the sign in front of the gluonic pole correlator. 

In the analysis presented here, we have tacitly assumed that the Wilson line in combination with the quark fields 
could be combined into a correlator with a gauge link and we only considered the cases of the simplest TMDs 
(without Wilson loops). To see what is happening in general we must realize that the complications arise from 
differentiating the Wilson line and obtaining covariant derivatives (in the above expressions the iD") in $g. This 
result is then separated into two parts, the parts <&g, which has the same color structure as $ itself and a gluonic pole 
part of which the color structure is Qq{x, x) — Tr c [^s(x)G na [p\], where the argument [p] of G na (see Eq. 173]) just 
indicates its (spatial) connection to the coordinates in the quark-quark correlator $(p), similarly as has been used for 
the Wilson lines U[p]. In the entangled expressions for the cross sections we have to use 

LC 



id? T . . . u l ^[ Pl } 4>( P2 ) u l : i] \pi] . 

= ($g(pi) . . . *(pa) . . .) + ($ 8 (zi) • • • 7rG na [ Pl ] $(pa) . . .) . (80) 

After further color disentangling of the right-hand-side, where one must be careful because ^(p 2 ) and possible other 
entries can have an octet structure, one then can recombine parts into tt^>q(x,x) = Tr c [$ 8 (x) irG na [p]], as explained 
for a full quark-quark-gluon correlator following Eq. [21 For an antiquark correlator ^{p 2 ) one obtains 

w? T ($(pi) . . . f/I n]t bi] l(pa) c/I" 1 Dpi] • • •) \ lc 

= ($g( Pl )...$(p 2 )...) - (^ 8 (x 1 )...^(p 2 )nG na \p 1 ]...y (81) 
Starting with the cross section o\ in Eq. [50] we find that weighting with transverse momentum gives 

<p? T oi) ~ Trc^g^!)^ A(^)rx] Tr c [$(a; 2 )r*A(z 2 )r 2 ] 

+ Tr c [$ 8 (x!) TrG na [ Pl } r* A(zi) r x ] Tr c [$(a; 2 ) T* A(z 2 ) T 2 ] 

+ Tr c [$ 8 (x!) rj A(z x ) r x ] Tr c [$(x 2 ) 2ttG"> 1 ] T* A(z 2 ) T 2 ] . (82) 
The traces in the last term are zero and we are left with 

{Pit a i ) ~ $%{x x ) + ^{x u x x ))§{x 2 )± x A(z x )A{z 2 ). (83) 
The weighted cross section of Eq. [57] gives 

{p? T a 2 ) ~ Tr c [$g(x 1 )nA(z 1 )r 2 $(a: 2 )r*A(z 2 )r 1 ] 

+ Tr c [$ 8 (xx) 7rG" Q [ Pl ] r* A(*i) T 2 $(z 2 ) T* A(z 2 ) r x ] 

+ Tr c [$ 8 (xx) rj A(zx) T 2 $(x 2 ) 2nG na \ Pl ] T* A(z 2 ) rj . (84) 
All terms survive and we are left with 

<Pi> 2 > ~ ($%(x x )+3tt$%(x x ,x x ))${x 2 )£ 2 A(z x )A(z 2 ). (85) 
Combining both contributions in the same way as in Eq. 1591 one obtains the result [ill ] , 

(Pit CT ) ^ (| ± 1 ^ 1+(D)1 (^i) Si - !v^t^ 1+D| (^) S 2 ) A(z0 A(z 2 ) 

= $g(x 1 )<I>(x 2 )(§±is 1 - 1 ^ T E 2 ) A(zOA(z 2 ) 

s — * — ' 

s 

+ n^(x x ,x x ) $(z 2 ) (^±i E x - ^ E 2 ) A(«i) A(z 2 ), (86) 
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where the combination of hard squared amplitudes in the second term of Eq. [86] is referred to as the gluonic pole cross 
section. In the case of quark-quark scattering with distinguishable quarks (Ei = £2) one thus has Eqp = N %_^ Ti, 
and for iV c = 3 we get 

(p? T tr) ~ ($g(n)$(a: 2 )E + 7r$S(n,a;i)$(a:2)|s) A(«i)A(2a). (87) 

In this result the two pieces in the cross section experimentally can be distinguished because of their time-reversal 
behavior. For instance single spin asymmetries have opposite time reversal behavior as compared to spin-averaged 
or double spin asymmetries. This time reversal behavior affects the paramctrization of the gluonic pole matrix 
element (containing T-odd distribution functions such as Sivers and Boer-Mulders functions). These functions appear 
convoluted with the gluonic pole cross section rather than the standard partonic cross section. For hadron-hadron 
scattering the partonic cross sections and gluonic pole cross sections have been tabulated in Ref. (39j . 

At this point we want to comment on the usefulness of our diagrammatic approach starting with an assumed 
convolution of soft correlators <$>(x,p T ) (integrated over p ■ P) and hard amplitudes. For this we have to realize that 
the operators involved have particular canonical dimensions d and twist t. Assignment of definite twist is only possible 
for local matrix elements or collinear correlators, e.g. the quark-quark (or gluon-gluon) corr elators$[*= 2 l(a;), & t=3 \x), 

\t — 31 \t — 31 

... or quark-quark-gluon correlators $ D (x) and $> A (x) . We have schematically 

$[ d=2 l(x,p T ) =$^ 2 ](lE,p T ), 

indicating on the right-hand-side that operators are involved of arbitrary twist with the minimal twist being the 
canonical dimension. 

Now, let us look at an observable at a high-energy scale Q including its dependence on collinear fractions and 
transverse momenta assuming appropriate identification of parton variables with kinematical variables, which for 
transverse momenta in hadron-hadron scattering would be q T = p lT + P2 T - We have 

da(xi,x 2 ,q T ;Q) = $ [d=2] (xi,p lT ) ® <S> [d=2] (x2,P2 T ) ® E 
M 

+ tt^^i^it) ® $ [d=2l (* 2 , P2T ) ® E] 

+ i[$ [d=31 (zi,Plr) ® <S> [d=2] (x2,P2 T ) ® E + $l d=2 ](n,p lT ) ® <S> [d ^(x2,P2 T ) ® ±] 

+ (88) 

of which only the first term survives at Q — >• 00. Looking at the gr-averaged (i.e. integrated) cross section, we have 
as leading contribution 

(a)(x 1 ,x 2 ;Q) Q -^° $' t=2 '(xi) ® ^ t=2 \x 2 ) ® E, (89) 
which summarizes the collinear approach. The leading contribution in the g T -weighted cross section is given by 

(q T a)(x 1 ,x 2 ;Q) Q -^° (p lT $ [d=2l )(a; 1 ) ® $ [ * =2] (x 2 ) ® E + ^^(x,) ® {p 2T $ [d=2] )(x 2 ) ® E, (90) 
where the weighted TMD correlators contain actually twist three and higher twist operators, schematically 
(p T ^ d=2 ^){x) = $[^ 3 l(x). 

Among the twist three correlators we have T-even correlators $^~ 3 \ $^ -3 ' in the combination $g -3 ' and the T-odd 

gluonic pole correlator ^q -3 ' (x, x). The usefulness of the diagrammatic approach starting with Eq. [S5]is in its ability 
to provide us in a straightforward way with the tree-level coefficients of the relevant (combinations of) twist-3 collinear 
correlators through the intermediate step of TMD correlators with a complex (process-dependent) gauge link. Similar 
results for the weighted asymmetries can also be obtained using collinear functions including the gluonic pole or ETQS 
functions directly from the start [4(| [4l[ . The tree- level results of the diagrammatic approach, of course, need to be 
improved upon by including NLO QCD contributions. In line with our proof of absorbing all gauge connections into 
a gauge link in 1PU processes, the diagrammatic approach also provides the tree-level results for higher transverse 
moments as long as only one initial state hadron is involved. This is useful first of all for leptoproduction processes 
involving a single incoming hadron and, furthermore, for selected asymmetries in hadron-hadron scattering. 
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B. Double weighted asymmetries 



We want to use the analysis in terms of transverse moments to indicate what is happening when the transverse 
momentum of two hadrons is involved, e.g. a double Sivers or double Boer-Mulders effect. The weightings with PxtPxt 
or with j> 2t can straightforwardly be obtained from Eq. I59[ although the splitting of TMDs with more complicated 
gauge links in T-even and T-odd parts is certainly not trivial. We will investigate the weighting with pf T following 
the same steps as done in the above for the single moments. 

We will first turn to an 'elementary' process, namely DY scattering. Having one color loop, one finds that in single 
and double weighting one needs 

Tr c [7] =-^Tr c [I]Tr c [J], 



1 



Tr c [T a T a ] = — Tr c [T a T a ] Tr c [J], 



Tr c [T a T b T a T b ] = - 



N C (N? - 1) 



Tr c [T a T a ] Tr c [T T ], 



yielding 

/ a 



$g {x X ) $g {X2 ) " TT^ (X 1 , X! ) $g {X2 ) " $g (Xl ) TT^ G (X 2 , X 2 ) 
~ lvJ3T 7r$ G'( a; l' X l) 7r$ G( a; 2,a;2)^ &DY- 



(91) 



This indicates a breaking of universality for the double gluonic pole contribution. This breaking is consistent with 
non-factorizability of DY at twist four (actually double twist-3), which would be the level needed to consider similar 
asymmetries in a collinear treatment that employs gluonic pole matrix elements and where one ought to find the 
same breaking. It also shows that TMD factorization of DY does not hold for the T-odd parts in the quark-quark 
correlators. The actual outcome is a strong suppression of the double T-odd contributions and a sign change (-1/8 as 
compared to 1). 

In the same way as for DY, we can analyze higher moments in other processes for which we will consider the 
example of quark-quark scattering. From the unintegrated result in Eq. [49] one gets the double- weighted result, 

(p? T J&*i) ~ Tr c [$g(a; 1 )r i A(2 1 )r 1 ]Tr c [^( 2 : 2 )r;A(z 2 )r 2 ] 

i r; A(*i) r x ] Tr c [$ 8 (* 2 ) -nG n % 2 \ r* a(z 2 ) r 2 ] 

nG na \p x ] rj A(zi) Ti] Tr c m(x 2 ) r* A(z 2 ) T 2 ] 
2wG n ^\p 2 ] rj A(z x ) T x ] Tr c [$ 8 (x 2 ) r* A(z 2 ) T 2 ] 
T* A(z x ) r x ] Tr c [$%(x 2 ) 2nG na \p x ] T 2 A(z 2 ) T 2 ] 
\ {nG na [ Pl ], nG n0 [p 2 ] } rj A (z x ) T x ] Tr c [$ 8 (x 2 ) r 2 A(z 2 ) r 2 ] 
p 2 ] Mxi) irG na [pi] H A( Zl ) T x ] Tr c [$ 8 (x 2 ) T* 2 A(z 2 ) T 2 ] 

r; A( 21 ) ri] Tr c [$ 8 (* 2 ) \{^[ P2 UG na [ Px ]}v* 2 a(z 2 ) r 2 ] 
r; A( Zl ) r x ] Tr c [nd na [ Pl ] $ 8 (z 2 ) nG^[ P2 } r* a(z 2 ) r 2 ] 



Tr c 


[*e(*i 


Tr c 




Tr c 


n(xi 


Tr c 




Tr c 




Tr c 


7rG™ ,3 [ 


Tr c 


'$&(xi 


Tr c 




Tr c 




Tr c 




Tr c 


$8<>1 


Tr c 


nG nl3 [ 


Tr c 


\G nSi \ 



■KG na \p x \Y\A{z x )Y x ] Tr, 
vG n % 2 \TlA(z x )Y x \T* t 
^G nf} \ P2 \Y\A{z x )Y x \T Xt 
p 2 \^>%{x x )Y\A{z x )Y x \ Tr, 



,[<i> 8 (x 2 )^G"' 3 [p 2 ]r*A(z 2 )r 2 ] 

[$ 8 (x 2 )7rG nQ [ Pl ]r*A(z 2 )r 2 ] 

[7rG na [p 1 ]$ 8 ( a;2 )r;A(z 2 )r 2 ] 
[$ 8 (x 2 ) 7 rG nQ b 1 ]r;A(z 2 )r 2 ] 

[i,G na \p x \^{x 2 )Y\A{z 2 )Y 2 \. 



(92) 



\p 2 ] $&{x x ) T* A(zi) Ti] Tr c 

Several terms are trivially zero after color tracing (four through nine). The others give particular transverse moments, 
where in particular for the last terms one must evaluate the extra color factor because of the specific color flow. To 
be specific for terms eleven through fourteen one needs 

Tr c [T a T b ]Tr c [T b T a ] 1 



Tr c [T a T a ] Tr c [T b T b ] - 1 ' 
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as compared to term ten. The result then is 

(p? T pL°i) ~ ^(x 1 )^x 2 ) + ^(x 1 )n^ G (x 2 ,x 2 ) + 7r^ G (x 1 ,x 1 )^ d (x 2 ) 

+ ^Uxx, xx) ir<S> G (x 2 ,x 2 )) Sj A( Zl )A(z 2 ). (93) 
From the unintegrated result in Eq. [55] one gets 

(p? t pL°2) ~ Tr c [ig( a;i )nA(z 1 )r 2 ^(x 2 )r;A^ 2 )r 1 ] 

+ Tr c [$g(xO r; A(«i) r 2 $ 8 (x 2 ) ttG"' 3 [ P2 ] r* A(z 2 ) r x ] 
+ Tr c [$8^) TrG na [ Pl } rj A(zi) r 2 $^(x 2 ) A(z 2 ) r x ] 
+ Tr c 2^g"> 2 ] r; A(zO r 2 $ 8 (x 2 ) r; A(z 2 ) r x ] 

+ Tr c [$ 8 (ari) rj A(*i) Y 2 W d {x 2 ) 2nG na [ Pl ] V* A(z 2 ) rj 
+ Tr c [$ 8 (zi) \{i:G na [ Pl ], nG n0 \p 2 }} T\ A( Zl )T 2 $ 8 (x 2 ) r* A(z 2 ) r x ] 
+ Tr c [7rG^[p 2 ] $8 ) TtG na [px] rj A(zi) T 2 $ 8 (z 2 ) r 2 A(z 2 ) r x ] 
+ Tr c [$ 8 (. 1 )r;A(z 1 )r 2 $ 8 (a; 2 ) ^{7rG n/J |p a ] s 7rG na [pi]}r5 A(z 2 )rx] 

+ Tr c [$ 8 (zi) r; A(zx) r 2 ttG">i] $ 8 (x 2 ) 7tG^[ P2 ] r; A(z 2 ) r t ] 

+ Tr c [* 8 (xi) 7rG na [pi] rj A(*i) T a $ 8 (x 2 ) 7rG^[p 2 ] T* A(z 2 ) rj 
+ Tr c [$ 8 (zi) 7rG n/3 [p 2 ] r* A{ Zl )T 2 $ 8 (z 2 ) 7rG" Q [ Pl ] T* A(z 2 ) r x ] 
+ Tr c [$ 8 (zi) 7rG n/3 [p 2 ] rj A(z x ) T 2 ttG""^] $ 8 (x 2 ) T* A(z 2 ) r x ] 
+ Tr c [irG^lpi] $ s (xi) rj A(z x ) T 2 $ 8 (^ 2 ) ttG">i] T* A(z 2 ) r x ] 

+ Tr c [7rG^[p 2 ] ^(m) rj A(z x ) T 2 7rG">i] $ 8 (:e 2 ) T* A(z 2 ) I\] . (94) 

To rewrite everything in terms of standard transverse moments $g we have to compare the color structure of the 
terms with two gluonic poles with the standard color structure for a 2 . We need the ratio 

Ti c [T a T b T a T b ] _ 1 
Tr c [T a T a T b T b ] ~ ~ 7V C 2 - 1 ' 

We note that the factor Tr c [T a T a T b T b ]/Ti c [T a T a ] Tr c [T b T b ] = l/N c is already incorporated in the factors multiplying 
Si and E 2 . The result is 

(PitpL^) ~ {n(x 1 )^x 2 ) + 3n{xi)^(xa,x i ) + 3^(x 1 ,x 1 )^x 2 ) 

+ TTfcgOEi, Xl ) tt$§(x 2 , x 2 j) t 2 A( Zl ) A(z 2 ). (95) 

Combining the terms we get 

(Pit P 2 t a ) ~ *g(x 1 )$g(x 2 )(g±iS 1 - ] ^ =T X: 2 ) A( Zl )A(z 2 ) 

\ ✓ 

+ ($g(xx) 7r$&(a; 2 ,a: a ) + ir^faxt) *g(a: a )) E x - ^ E 2 ) Afo) A(z 2 ) 

V v ' 

+ Tri&gOci, Xl ) n$ G (x 2 ,x 2 )( ( N \%f_% +3) Ex - i^Izi| s 2 ) A(z x ) A(z 2 ), (96) 
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which in our example of distinguishable quarks (Si = £2) gives 

(Pit P2t a ) ~ ^(x 1 )^ d (x 2 )tA(z 1 )A(z 2 ) 

+ {n(x 1 )7T^ G (x 2 ,x 2 ) + 7r#g(ii, Xl ) $^(x 2 )) *§± t A( Zl ) A(z 2 ) 

+ 7r$g(zi, zi) 7T$S(^, x 2 ) (N \-fi^ 21) t A(zx) A(z 2 ) (97) 
= (%0n) $|(a: 2 ) E + (4g(xi) 7r$^(a; 2j i2) + ^(n.n) ^(z 2 )) \ t 

+ ir$%{x x ,xi) ir<5> G {x 2 ,x 2 ) i| s) A( Zl ) A(z 2 ). (98) 

The result shows that while the effect of a single gluonic pole is suppressed (factor 1/2) as compared to the naive 
TMD factorization, the effect of the double gluonic pole (double Sivers or double Boer-Mulders) is suppressed less 
than the naive expectation expected on the basis of single gluonic pole effects that are expected to show up in single 
spin asymmetries (factor 15/32 versus 1/4). 

Although these examples for DY and hadron-hadron scattering show that TMD factorization fails, we think that 
the starting point Eq. 25] at the end of Section IIIII appears to be useful, incorporating all leading matrix elements 
based on a counting of canonical dimensions. From that expression, one can proceed and calculate higher weighted 
asymmetries. 



VI. CONCLUSIONS 



In this paper we have used a diagrammatic approach to analyse the leading contributions at tree-level in hard 
processes in which several hadrons are involved. The diagrammatic approach combines correlators involving hadron 
states and parton fields with hard partonic amplitudes. The correlators depend on (on-shell) hadron momenta (initial 
state hadrons P, produced hadrons Kh) and the parton momenta (p and fc). At high energies, momentum fractions 
1? or (l/z)Kh and transverse momenta can be used, while the components p-P and k-Kh are integrated over. 
At that stage one has correlators that also have a natural link to light-cone wave functions [42|, [43| . In the usual 
collinear treatment, the transverse momenta are also integrated over, or one studies the limit in which these momenta 
become large, which just is a collinear treatment involving a partonic process with one or more additional partons. 
Using TMD correlators one aims at incorporating all features related to (soft) transverse momenta, including in the 
parametrization of TMDs the possibility of T-odd correlators, not forbidden by any symmetry. These incorporate the 
effects of initial and final state interactions [44| and they can explain the existence of single spin asymmetries at high 
energies. 

The basic tree-level result was presented in Section IIIII It includes all Wilson lines originating from collinear 
gluons in all of the correlators. These contributions are sufficient to study the collinear correlators (integrating over 
transverse components of the parton momenta) and give rise to straight light-like gauge links. For TMD correlators, 
one has to include transverse gauge connections, but after accounting for the collinear gluons, this all works out nice 
and produces the transverse pieces at the right place [8, 9]. One obtains a color gauge-invariant expression (Eq. |4"2"|) . 
which is fully entangled. 

In order to disentangle color, one needs to consider for each diagrammatic contribution in the hard amplitude, 
also its color-flow possibilities. This does not resolve the entanglement in general. Also, integrating over transverse 
momenta of final state partons, e.g. when one considers jet production, does not resolve the problem. Then one 
remains with an entangled situation if the initial state involves more than one hadron, as in the case of hadron-hadron 
scattering. This was considered as an explicit example of factorization breaking in Ref. [3?} ■ The full entangled 
tree-level result for this situation is given in Eq. 221 Even in cases in which one just has a single color flow such 
as in the electroweak processes of leptoproduction, electron-positron annihilation or the Drell-Yan process, one gets 
TMD factorization with correlators having specific process-dependent gauge links running via ± light-cone infinity. 
For these cases one then has a starting point with universal correlators and one can study the full QCD factorization. 
There are several examples of similarly simple processes, such as two-photon production in Drell-Yan, which even for 
gluons has a simple color flow (45|. 

The main aim of this paper was to show what happens in the special case of a 1PU process, in which only the 
transverse momentum in one specific hadron is considered. In that situation, one can resolve the entanglement (at 
least at tree-level). What remains is a color gauge-invariant expression (Eq. 159")) with a process-dependent gauge 
link. All-but-one collinear processes may seem hard to realize at first sight, since the parton transverse momentum 
is an integration variable, rather than an observable. Nevertheless, symmetry considerations in combination with 
polarization, in particular transverse polarization of target hadrons, may help to create the conditions for a 1PU 
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treatment. Furthermore one has in leptoproduction processes a whole class of such processes. Consideration of TMD 
correlators is possible for quark as well as gluon correlators. A recent example of the latter situation was pointed out 
in Ref. (4(| which considered heavy quark production in the Drell-Yan process. 

The simplifications and steps towards a TMD factorizable result become more transparant if one constructs trans- 
verse moments. This analysis has been given in Section[V]for single and double weighted asymmetries. Single weighted 
asymmetries can use the results of the 1PU situation and lead to the results where T-even transverse moments are 
convoluted with the standard partonic cross section and T-odd transverse moments are convoluted with a different 
(but color gauge-invariant) combination of squared amplitudes, the so-called gluonic pole cross sections. The T- 
odd transverse moments are weighted TMD correlators, involving p^,-weighted T-odd distribution functions like the 
Sivers or Boer-Mulders distribution functions of which the operator structure is a gluonic pole matrix element, i.e. a 
quark-quark-gluon matrix element at vanishing gluon momentum. 

While the single weighted asymmetries exhibit factorization at the TMD-level (with process-dependence through the 
link in the TMD) and at the collinear level involve gluonic pole cross sections, rather than the standard partonic cross 
sections, the explicit evaluation of the double weighted asymmetry shows that such factorizations are not generally 
applicable. Additional color factors come in at this level, which provide besides double weighted TMDs for hadron 
1 and hadron 2, products of single weighted functions. The factors multiplying the final result for the cross terms, 
however, is not related to the factors appearing in the single weighted results. It does involve products of the T-even 
and T-odd transverse moments, but each of these terms with their own factors. 

The analysis in terms of transverse moments, however, also provides us with an important result for correlators 
belonging to final state hadrons. In analogy to the analysis of distribution correlators, the transverse moments for 
fragmentation correlators involve a part Ag with a T-even combination of operators and a part Ag with a T-odd 
combination of operators, the latter involving an integration over transverse gluon fields corresponding to vanishing 
gluon momentum. Such gluonic pole matrix elements (and also multi-gluonic pole matrix elements appearing in 
higher transverse moments) vanish for fragmentation. This has been shown in models as well using spectral analysis 
of final states [13, HH or studies of analytic properties [29] based on the field theoretical correspondence to amplitudes 
mentioned at the end of section [TTJ Since the color structure of the remaining Ag part is simple, there are no 
complications to absorb the appropriate gluons into gauge links. This implies that the result of Eq. also applies 
for unintegrated situations, where the dependence on transverse momenta in the final state is kept. Thus we can use 
in Eq.[49]TMD correlators A(z,k T ) rather than the collinear correlators A(z). Instead of Eq.|42]one has a tree-level 
result of the form 

da ~ Tr c [[/l[p 2 ] $( Xl ,p lT ) U-[p 2 ] Wl[ Pl ,p 2 ]]Tr c [ul[ Pl ] <5>{x 2 ,p 2T ) U.[p{\ Wl[p u p 2 }} 

x r* A( Zl , k lT ) r 1 r* A(z 2 , k 2T ) r 2 , (99) 

with entanglement only involving the initial state TMD correlators, including dependence on color flow in initial and 
final state. In spite of having a T-even operator structure, the universal TMD fragmentation correlators A(z, k T ) 
or its transverse moment Ag(z) allows for T-odd fragmentation functions in its parametrization such as the Collins 
function because the non-plane- wave nature of the states \Kh,Kx) in the matrix elements of Eq. [5] prevents using 
time-reversal symmetry constraints. 

The formalism in this paper allows a rich phenomenology, including many results on single transverse momentum 
weighted cross sections that have already been obtained. It can be used for quark as well as gluon TMDs, although 
the latter involve a number of additional complications. However, we do want to emphasize once more that our results 
are strictly tree-level, which does provide insight into the universality of TMDs and which is a necessary condition to 
study factorization issues, but it does not provide proof of factorization in QCD. For a discussion of those aspects we 
refer to Ref. (T^j, to Ref. [l6| on subtleties with gauge link structures at higher orders and to Ref. [l8| for attempts 
to bridge the gap between the phenomenology and the more formal definitions. 
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Appendix A: Collinear gauge links 

In this appendix we incorporate the collinear gluons, i.e. the first term in the Sudakov expansion of A fJ, (p) given in 
Eq. [2U This will produce the parts of the gauge link along light-like directions conjugate to the momentum P and 
in the collinear situation (when all transverse momenta are integrated over) it produces the gauge link as given in 
Eq. [17] in the light-cone correlators. 

Looking at the Fourier transformed fields with collinear Wilson lines starting at ±00, we consider the Fourier 
transform of the field ^ For the field with Wilson line starting at minus infinity we get, 



E/L n V(p) 



d 4 £ exp (ip- i) exp 



-w 



d(r] ■ P) n- A(rj) ip(x). 



(Al) 



Looking only at the relevant component p + — p ■ n one gets, 



u^Ap) = £(-»; 



N 



N=0 



di] 



N 



dr\ 



N-l 



drh A+(^)...A+(»jrMOe ip+r 



(A2) 



where the arguments run between —00 < r/ N < rj N _ 1 < . . . < r q 1 < £ , implemented through 6 functions 6(rj N1 
r]]j) . . . 6{rf[ — ?7^)#(£~ — r li)^ which can be rewritten as momentum-space integrations, 
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Including again the dependence on other momenta, we find 

d A p N f d A Pl A n (p N ) A n ( PN ^- PN ) 
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(2tt)4 
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(2tt)< 
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tp(p-pi) 



d A PN A n ( Pl ) 

(27r) 4 (xi + ie) (xi + X2 + ie) (xi + . . . + xn + ie) 




where Xj — Pi ■ n. We thus find for the gauge link 



M=0 



(A3) 



where = 1, The gauge link and the terms in its expansion not only have a particular structure in coordinate or 

momentum space, but they also have a charge structure. In particular for applications in non-abelian gauge theories 
one has matrix- valued fields A^ — A^ a T a . 

The gauge link can be written in a nicer symmetric form for the correlators. Commutators [A n (pi) , A n (P2)] don't 
matter here, neither in color space, where they are 'contained' in other matrix elements with less gluons, nor in Hilbert 
space, where one has commuting fields, which for fields with the light-cone index n imply vanishing commutators. We 
can then use relations 



(xi + x 2 + ie) 



1 



1 



(xi + ie) (x 2 + ie) 
1 



1 



(xi + x 2 + x 3 + ie) I (xi + x 2 + ie) (xi + ie) 



(xi + ie)(x 2 + ie) ' 
1 

(xi + x 2 + X3 + ie) 



(A4) 



(xi + ie)(x2 + ie) 



6 permutations 



3 permutations 



(xi + ie)(x 2 + ie)(x 3 + ie) 



(A5) 
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k-Pi-p' 




k-pi k 



k-Pi-p'i 




P-Pi 



k-Pi-p'i 



(a) 




P-f ^ 
(b) 




p-pi 



P1+P1 



(c) 



FIG. 9: The gluon insertions on an outgoing quark line coming from two different soft pieces, one from $(p) and one from 
respectively. 



and its generalization to more terms, to symmetrize the result. This gives 
,[p)= I <i 4 £ exp (i v ■ £) & exp ( -iq I 



U [ l l] ijj(p) = / d 4 £ exp (ip ■ £) ^exp ( -ig f d{ri ■ P) n ■ A{rj) I tp(x) 



sr±[^PL [ d4 P N An (pi ) A " M- A " (p^y ) / f y ^ (am 

Nl J (2tt) 4 ' ' ■ 7 (2tt)4 (xx + ie)(a: 2 + ie) . . . (x N + ie) V [ P ^ Pl j ' 1 ° j 



AT=0 



From this expression one sees that the term U^^ M ^ is the consecutive action of M simple (commuting) U± - 
connections 

For a link along n coming from +oo one has 



j4. 



LTf^O) = J d 4 £ exp(ip-£) ^exp^-i 5 J d(rj ■ P) n ■ A(rj) J tp(x) 

(2tt)4"7 (2^) 4 (-^ + ie) (-s w _i+ie) ^ (-Xi+ie) W Plj 



(27r)4"7 (2ir)*(-x 1 +ie)-(-x 1 -...-x N +ie) iP \ P j£\ 



fi/^i / ^ ^"(Pi) ^"(pa) ■ ■ ■ A n { PN ) ( A \ 

^ TV! 7 (2^) 4 ' • 7 (2tt) 4 (-n + ie)(-x 2 + ie) • . . (~x N + ie) V \ P ^ Pl \ ' 



_ 1 f d 4 Pl f d 4 PN A n (p x ) A n ( P2 ) . . . A n ( PN ) 

\ ■ n 



The notation used in equations with many fields and links will be U + ip(p) — U+ [p]V"(p)- 



Appendix B: Intertwined gauge connections 

We consider the situation of gluon insertions on an outgoing quark line coming from two different soft pieces, one 
from <&{p) and one from $(p'), respectively. There are three leading contributions of A k (pi)p^ and A k (jp' l )p-y gluon 
components, corresponding to the diagrams in Fig. [51 one of them involving a three-gluon vertex. The combined 
result of the diagrams actually nicely adds up to 

A n = \${k)U^%,p>]Ti,{pj\ ...^{p') =^(k){U^%'],U^\p}}T^(p)] ...^{p'l (Bl) 

which is a gauge connection which is the (color) symmetrized product of simple connections. 

In order to illustrate the recursive procedure, it is instructive to give the result for the situation of insertions on 
an outgoing quark line for three gluons coming from two different soft pieces, two gluons coming from $(p) and one 
coming from $(p'), respectively. This is shown in Fig. [10] and involves three- and four-gluon couplings. The result 
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P-P1-P2 




(e) 




(f) 





P-P1-P2 (b^/pi rt) u/ Pi 
(c) 




FIG. 10: The gluon insertions on an outgoing quark line for three gluons 
coming from two different soft pieces, two gluons coming from &(p) and 
one coming from respectively. Note the absence of a three-gluon 

vertex with gluon legs with momenta p\ and pi coupling to the same 
soft part $(p). Note the absence of a three-gluon vertex with gluon 
legs with momenta pi and P2 coupling to the same soft part <&(p). It is 
actually already included in the soft part with one gluon leg less. 



becomes 



[mu l * m) {p,p'}r...^(p)]...^(pi) 

mm ^ m \p] u^%'] + \ u^%] u^%'] u^%] 

+ \U^ {1) [p']U^ {2) [p\)T...^{p)\..^(p') 



\uf {1) [p]uf (1) [p']uf (1) [p] 



. ^KD [p/] ^KD [p] v im) M)r _ _ _ m] _ _ _ W) 



(B2) 
(B3) 
(B4) 



For a general A^fe-gluon term, with i-, j- and k gluons from three (or more) different correlators, one finds an 
expression where the A- fields from these different correlators are a color symmetrized product of f/WC 1 ) factors 
containing commuting gauge fields. Thus if one, as is the case in Eq. 1291 looks at the result from insertions on leg 
&2 coming from correlators 4>(j>i), $(^2) and A(fci) it can be broken apart in a symmetrized product of simple gauge 
links, 



U [ ^ ] [P1,P2M] 



S{U^ ] [ Pl }U^ ] [ P2 ]U^ 



[ki]h 



(B5) 



of which the ordering is irrelevant. 
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